Shape and symmetry of a fluid-supported elastic sheet 
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A connection between the dynamics of a sine-Gordon chain and a certain static membrane folding 
f^f-^ . problem was recently found. The one-dimensional membrane profile is a cross-section of the position- 

time sine-Gordon amplitude profile. Here we show that when one system is embedded in a higher- 
(^ I dimensional system in this way, obvious symmetries in the larger system can lead to nontrivial 

^SJ ■ symmetries in the embedded system. In particular, a thin buckled membrane on a fluid substrate 

, , ' has a continuous degeneracy that interpolates between a symmetric and an antisymmetric fold. We 

O ,' find the Hamiltonian generator of this symmetry and the corresponding conserved momentum by 

_^^ , interpreting the simple translational symmetries of the sine-Gordon chain in terms of the embedded 

^ ' coordinates. We discuss possible extensions to other embedded dynamical systems. 
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PACS numbers: 46.32.-|-x 46.70.-p 68.60.Bs 81.16.Rf 89.75.Kd 



I. INTRODUCTION 



'"O . Thin elastic sheets form rich patterns when stressed or confined. One may try to account for these patterns as 

G ■ weak deformations of a simple (e.g., flat) reference state. It has been recognized, however, that as the sheet is made 
increasingly thin and bendable, the validity range of such near-threshold analyses shrinks indefinitely, along with their 
experimental relevance jl|-|3| . Crumpled paper is a prime example of such a far-from-threshold behavior |4J] . 

A thin elastic sheet, compressed on top of a heavy fluid |5l-[l5|. provides a simpler example of this situation. The 
leading-order analysis of the deformation yields a periodic pattern of wrinkles with wavelength A = 27r[B/{pg)]^^^ 
[a, Q . The finite wavelength arises from a competition between the resistance of the sheet to bending (B being the 
r-o ■ sheet's bending rigidity) and the resistance of the fluid substrate to height changes (p being the fluid mass density 
^\ ' and g the gravitational acceleration) . Since B depends on the sheet thickness w as B '^ w^ , the wrinkling wavelength 
decreases with decreasing thickness as A ~ w^'"^. A higher-order analysis reveals that the wrinkling pattern is actually 
stable only for small lateral displacements, A < Aw ^^ A^/i, where L is the total length of the sheet [ISl, [ij|. This 
f~^ , validity range diminishes with either decreasing w or increasing L. Beyond it the deformation becomes localized 
Cf^ ■ in a finite domain (fold) of width k~^ ^ A^/A [ij, |lj|. It takes a mere displacement of A ~ A to concentrate all 
the deformation into a single wavelength and make the fold contact itself [^, [ij, [l^. Thus, the relevant response 
of a sufficiently thin sheet is practically always far from threshold, involvin g a strong localized deformation. Similar 
observations hold when the fluid substrate is replaced by an elastic one [iJ: KL6l - l22| . 

It is therefore fortunate and beneficial that the fluid-supported sheet under uniaxial compression turns out to be 
an integrable system [15]. Exactly solvable problems often owe their integrability to extra symmetries. In the current 
work we focus on the corresponding extra symmetry in the fluid-supported sheet, which is manifest in a continuous 
degeneracy of ground states |23l |2Jj . We connect this peculiar property of the one-dimensional buckling problem to 
the simple translational symmetries of the two-dimensional sine-Gordon problem. 

In Sees. lU and Hill we repeat the deflnition of the problem and its representation as a dynamical system [l^. In 
Sec. lIVI we establish the relation between our problem and the sine-Gordon chain, and in Sec. |V]we characterize the 
resulting extra symmetry. Section IVII describes consequences concerning the fluid mass displaced by the sheet and 
the associated buoyant force. Finally, in Sec. lVIII we discuss consequences of the new symmetry and analogous cases 
in other physical systems. 
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II. THE SYSTEM 

The system under consideration is schematically illustrated in Fig. [T] An incompressible sheet of bending modulus 
B lies on a semi-infinite fluid of mass density p. The sheet is assumed to be indefinitely long and deform in the xz 
plane, remaining uncurved along the y axis. The deformation is parameterized by either the height profile h(u) or the 
slope angular profile 4>{u) as a function of the arclength u along the sheet. The two functions are geometrically linked 
by the relation h = sine/), where hereafter a dot denotes a derivative with respect to u. The sheet is subjected at its 
edges to a uniaxial pressure P along x, which is accompanied by a displacement A = J_ du{l — cos </)) between the 
two edges. 




FIG. 1. Schematic view of the system. 



A given deformation costs energy per unit length, E — Eb + Eq, with contributions from bending, Eb — 
{B/2)J_ du(p^, and from gravity, Eq = {K/2) J_ duh^ cos (f), where K = pg. We use dimensionless variables, 

where all lengths are scaled by {B/KY^'^ = A/(27r) and all energies by B. The pressure is consequently scaled by 
(BK)^'^. An equilibrium shape of the sheet is one that minimizes E~PA for a given P and the appropriate boundary 
conditions. We thus look for the minimum of the following functional: 



S 



C 



1 



12 



duC{4>, h, (j), h) 
1 



+ 



:/i2, 



P(l 



Q{su\(f) — h). 



2" 2 
In Eq. ([l} another Lagrange multiplier, Q(u), has replaced the local constraint on the relation between h and 
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III. THE SHEET AS A DYNAMICAL SYSTEM 



We study the deformation of the sheet using a mechanical analogy, where the arclength u along the sheet stands 
for time, and the angle 0(m) and height h(u) represent two coordinates with conjugate momenta p^ and ph Il5| . 
The equilibrium shape of the sheet is given by the trajectory [(/)(m), h{u)] that minimizes the action in Eq. ([l}. The 
Lagrangian £(</>, h, 0, h) appearing in Eq. ([Ij transforms into the following Hamiltonian: 
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n{(t),h,p^,Ph) = -pl+phsm(f>- -h'^cos(p + P(l 
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The fact that H is a constant of motion corresponds to one apparent symmetry of the homogeneous sheet — its 
invariance to translation in u. In what follows we explicitly distinguish between derivatives with respect to the 
canonical variables (using the symbol d) and derivatives with respect to the coordinates (using dots and primes). 
Hamilton's equations yield the following dynamical system: 



h = dp,^ % = sin 4> 
P4, = -d^T-L = -ph cos - (ft.^/2 + P) sin 
Ph = -dhU = h cos 0, 



(3a) 
(3b) 
(3c) 
(3d) 



whose solution gives the equilibrium shape. Equations (|3a)) and (j3dp identify the ^-conjugate momentum as the local 
curvature, p^ = 4>, and the ^.-conjugate momentum as the fluid mass displaced by the deformed sheet up to point 
u, ph = /_ duihcoscf). In the case of a localized fold in an infinite sheet the following boundary conditions apply: 

/i, 0, = P0 = at u — >■ ±cxD, resulting mH = 0. Further requiring that = far away from the fold implies, through 
Eq. ([He]) , that also Ph = a.t u —^ ±oo. We have found the exact solution of this problem by relating it to the integrable 
sine-Gordon chain [ISj . The integrability of our system implies the existence of a complete set of conserved canonical 
momenta. Since the system has two canonical momenta, it must thus have two conserved dynamical variables, the 
Hamiltonian H and one other. We shall call this conserved quantity JC and determine it in Sec. |Vl We need to discuss 
first the relation between our system and the sine-Gordon chain. 

IV. RELATION TO THE SINE-GORDON CHAIN 

Equations dSa]), ^, and @ with H = 0, yield 

'^ + {<j)^/2 + P)^ + h = 0, (4) 

which is the known equation of Euler's elastica [4il25{. with the hydrostatic pressure term h {pgh in dimensional terms) 
playing here the role of an external normal force. Differentiating once with respect to u, we obtain a fourth-order 
equation for the angle (t){u) alone, 

"<j)' + [{3/2)<j)^ + P]4> + sm(j) = 0. (5) 

Equation ^ is related to a known hierarchy of nonlinear partial differential equations — the combined sine-Gordon 
and modified KdV (sG-mKdV) hierarchy [26|. The first three equations in that hierarchy, for the two-dimensional 
function (j){u,v), read: 

(/)' + 2(/3e''^ - ae-"^) = (6a) 

^' - i4> + 2{Pe"l' - ae-"^) = (6b) 

0' + {i/4)"(f>' + {i/8)4>^4> - ic4> + 2{l3e"^ - ae'"^) = 0, (6c) 

where a prime denotes a derivative with respect to v. Note that all equations in the hierarchy are invariant to 
translations in both coordinates, u and v. Setting a = /S = i/4 in Eq. (|6ap yields 

<i)'^sm(t). (7) 

This is the sine-Gordon (SG) equation in light-cone coordinates, u — {x + t)/2,v — {x — t)/2, describing the swaying 



angle 4>{x, t) of a pendulum at position x and time t along a chain of coupled pendulums. Setting in Eq. (I6bp 
a = /? = — 1/4 and projecting the equation onto a constant w, we obtain the physical-pendulum (PP) equation, 

-I- sin = 0, (8) 

where here 4'{u) describes the swaying angle of a single pendulum at time u. Finally, upon setting a = /3 = 1 and 
c — — P, and projecting onto a constant v, Eq. ([6cl) coincides with the shape equation for the sheet, Eq. ^. 

The relation to the sG-mKdV hierarchy enables us to obtain solutions of Eq. ([SJ from known solutions of the PP 
equation ([5]) or projected solutions of the SG equation ([7]). Specifically, we take the "breather" solution of Eq. ([7]), 



4>{u, v) = 4 tan 



K, sin(fc(M — w + ci)) 



fc^ + K^ = 1, (9) 



k cosh(K(u + V + C2)) 

describing a standing localized wave in the chain of pendulums. The two arbitrary constants, ci and C2, reflect the 
two symmetries of the SG chain under translations in space x and time t (equivalently, along the two light-cone 
coordinates u and v). Projection of the breather solution onto one of the light-cone coordinates (e.g., u with v — 0) 
yields the solution to the localized-fold profile, 

T \ K, sin(fc(u -|- Ci)) 1 o T , , 

0(w) = 4tan-i - :.\ J'' , k^ + K' = l. (10) 

\_k COSh(K(M -I- C2))\ 

Substituting it in Eq. ([5]), we find 

k^\{2 + Pf'\ n=\{2-Pf'\ (11) 

The expression for the height profile, h{u), can be obtained in closed form as well but is postponed until the end of 
Sec.lVl 
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FIG. 2. (color online). Shapes of a folded sheet for P = 1. In between the symmetric fold (lowest black curve) and the 
antisymmetric one (uppermost red curve) there is a continuum of degenerate shapes having the exact same displacement and 
energy. The shapes were calculated from Eqs. (|10[) and (|lip with C2 — and Cifc = 0, n/8, 7r/4, Stt/S, and 7r/2. 



THE EXTRA SYMMETRY 



The discussion above shows that the folding profiles 4'{u) are embedded in an autonomous two-dimensional system: 
the sine-Gordon system. The obvious translational invariance of the sine-Gordon system in v entails a corresponding 
invariance of the folding profiles that is not obvious: not only can the entire profile be translated along the sheet, but 
the oscillatory and decaying parts of the profile can be shifted independently [cf. Eq. (|T0)) ]. Consequently, in between 
the symmetric fold (ci — C2 = 0) and the antisymmetric one {cik = tt/2, C2 — 0) there exists a continuous family 
of degenerate profiles, having the exact same displacement and the exact same (minimum) energy |23l. |24|. This is 
demonstrated in Fig. [2j 



Conservation law 



To find the conserved quantity associated with the extra symmetry, lC{(j),h,p^,pfi), we seek the generator of v- 
translations in the Hamiltonian phase space of Eq. ([2|) , 



{/,/C} = /', 



(12) 



where f{u,v) is an arbitrary function expressed in terms of our canonical variables [Eq. ([2])], and {■ • ■ } are Poisson 
brackets. In the derivation below we make use of the following: (a) our knowledge of the first conserved quantity, H 
[Eq. ^], and the resulting Hamilton equations ([3|); (b) the set of equations, analogous to Hamilton's, that derive from 
K., i.e., {(/), /C} — dp^lC, {h,IC} = dp^lC, {p^,IC} — —d^lC, {p/i,/C} = —dhIC; (c) the fact that 4>{u,v) should satisfy 
the SG equation ([T]). (Note, however, that we make no use of any particular solution of the SG equation.) Below we 
obtain a set of differential constraints on /C, find a functional form consistent with these, and finally verify that this 
/C, (a) generates the desired shift in v, and (b) commutes with %. 

Taking f — (j) in Eq. (J12p and differentiating both sides with respect to w, we find {d/du){(j), 1C\ = 0' = sini^ = h. 
Thus, up to an integration constant, taken as zero. 



dp^JC = 4>' = h. 
We note that the equation (j)' ~ h provides the connection between the coordinate v and the actual sheet, 



(13) 



hdv = (j)du. 

It also represents the fact (now in terms of v rather than u) that (p and h are inter-independent; just as h = sin (j) along 
the u coordinate, (/)' — h along the v coordinate. Applying a similar procedure with f = h yields {d/du){h, IC} = h' = 
{siiKpy = cos <p(j)' = hcoscj) = ph- Thus, up to an integration constant. 



dp,)C 
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Next, for f — p^, we have {p^, /C} — p'. — ip' ~ sin </>. Thus, 

d^lC = — sin0. (15) 

Finally, for f = ph, we get {d/du){ph, K,} = p'y^ = [h cos 0)' = cos 4)h' ~ h sin 00' = ph cos — /i^ sin 0. Using Hamilton's 
equation ([5cl) . we write this result as ((i/d?i){p/i,/C} = —p^ — [(3/2)/i^ + P] sin0 = —{d/du)[p^ + h^/2 + Ph]. Thus, 
up to an integration constant, 

dhlC^P4, + ^h'' + Ph. (16) 

Gathering the information from Eqs. (|13 p - (|16p . we obtain an expression for the second conserved quantity, 

IC = \pI + hp^ + ^h^h^A + P) - (1 - COS0). (17) 

We may now verify that this candidate form for /C satisfies the requirements (a) and (b) above. By construction, 
the K. of Eq. (IT7)) satisfies Eqs. p^ ~([TC l) . Thus, it generates a shift in v, fulfilling requirement (a). To address 
requirement (b) we may explicitly verify that /C of Eq. (|17p is a constant of motion, 

{IC,H} = {d^lC){dp^H) - {dp^lC){d^H) + {dhlC){dp,H) - {dp,IC){dhH) = 0, (18) 

as can be verified by direct substitution using Eqs. (|13p - P^ . While this /C bears no obvious resemblance to H of Eq. 
([2]), the two must nevertheless be equivalent in generating the fold. This is because a localized SG solution along a 
line of constant u must be equivalent to that along a line of constant v. To verify this duality we form the Hamilton 
equations generated by /C, 

(j)' ^ dp^JC :^ h (19a) 

h' = dp,IC^Ph (19b) 

P'4, = -d^lC = sin (j) (19c) 

p',, = -dhJC = -p^ - /1V2 - Ph. (19d) 

When the two dynamical systems ([3]) and (I19p are compared, their equivalence is not obvious. The replacement 
/i ■<-)• P0, however, makes Eqs. (J19al) and (|19cl) the same as Eqs. (|3a|) and (l3b)) . respectively, and Eqs. (J19bp and (I19dp 
yield Eq. ^ of Euler's elastica. Once a single equation for 0(w) is constructed, 

0"" + [(3/2)(0')' + ^]0" + sin0 = 0, (20) 

and compared with Eq. ([5]), the equivalence is apparent. 

Finally, we use the extra symmetry to obtain closed-form expressions which could not be derived in earlier works. 
The height profile, h{u) = J_ duismcj), is obtained according to Eq. (|19ap by differentiating the breather, Eq. ^, 
with respect to v and then setting v = 0. We get 



k cos ui cosh U2 + K sin ui sinh U2 
k"^ cosh U2 + K^ sin ui 

Ul^k{u + Cl), U2^k{u + C2). 



h{u) = -^AkK '' "" ,';^ ""7;:^ ^ " ^^^: "^^ :"" "^ (21) 



Similarly, according to Eq. (J19bp . we obtain the profile of displaced fiuid mass, pt — J_ duih cos cf), by differentiating 
the breather twice with respect to v and setting v = 0. This yields 

2kK 

Ph{u) = 2 2 \2kn{k^ cosh(2u2) + k^ cos(2ui) + fc^ — k^} cosui sinhM2 

(fc^ cosh 52 + K^ sin ui)2 

-{{k^ -K^)[fc2cosh(2u2) -I-k2cos(2ui)] + fc** + 6fc^K^ + k"*} sin ui cosh M2] . (22) 

B. Shape transformation 

The actual transformation generated by the extra symmetry is nontrivial. From Eqs. (|19ap and (J19bp we see that a 
slight translation 5v in the extra dimension adds to the angle of the sheet at each point a small amount proportional 
to the height at that point, 4>{u) -^ 4>{u) + h(u)Sv] to the height it adds a small amount proportional to ph, the fluid 



mass displaced up to that point, h{u) — )■ h{u) +pji{u)Sv. In Fig. [2]one can try to follow these small changes between 
consecutive curves. 

We verify that such an infinitesimal transformation leaves the mechanical energy of the sheet unchanged. Applying 
it to the bending energy, we get 

SEb 

— - — :^ (j)h = sin (fxj) = —(d/du) (cos (j)), (23) 

ov 

which is a total differential, contributing only a boundary term to the integrated bending energy. That contribution 
vanishes for our boundary conditions. For the gravitational energy we get, with the help of Eq. (I3dp . 

^ = hph COS0 - ift3 gin0 = phph - h^h = {d/du){pl/2 - hyS), (24) 

ov 2 2 

which is a total differential as well. The contribution from the /i^ term vanishes for our boundary conditions; the 
vanishing of the second boundary contribution is discussed in the next section. Thus, both the bending energy and the 
substrate energy are separately invariant to the transformation. We further confirm that the displacement is invariant 
too, 

^ = hsm(j} = hh = {d/du){h^/2). (25) 

ov 

Finally, the variations Sh/6v = ph and 5sui(t>/5v = hcoscf) are equal thanks to Eq. (jSdl) . ensuring that the transfor- 
mation does not violate the geometrical constraint h = sin (f>. 

VI. BUOYANT FORCE 

When defining the system in Sec. |lll we have prescribed the lateral force per unit length acting on the sheet, P, 
which is coupled to the lateral displacement, A. By contrast, no explicit term has been introduced to impose an 
external vertical force or a certain vertical displacement, apart from the boundary condition h{u -^ ±(X)) = 0. This 
is equivalent to imposing a zero vertical force. In the absence of such an external force, the total vertically displaced 
fluid mass must vanish. The localized solutions of the profile equation should automatically satisfy this constraint. 

As mentioned above, the mass of vertically displaced fluid up to a certain point u in the sheet is equal to the 
/i-conjugate momentum, ph = J^ duihcoscf) [cf. Eq. (j3dl) ]. Thus, the total displaced fluid mass is Mf — ph{oo). The 
fact that A/f = then readily follows from Eq. ([5c)) applied at u — >■ oo. It is inferred also from the conservation of 
/C when setting /C(— oo) = /C(oo) in Eq. (|17[) . Finally, we explicitly confirm that p/i(±oo) — ;■ using the closed-form 
expression for ph, Eq. (f^ . 

VII. DISCUSSION 

The physical property of the sheet-fluid composite, underlying the extra symmetry, is unclear. We have found no 
physical argument to explain why the shape transformation described in Sec. IV Bl leaves the mechanical energy of 
the sheet invariant. Nonetheless, in the two-dimensional problem, in which we have embedded our one-dimensional 
system, the extra symmetry is a trivial translation invariance along the extra dimension. This situation is reminiscent 
of the phason freedom in quasicrystals [27H29| . In the actual quasicrystal the degenerate phason states correspond 
to nontrivial coordinated flips of atom locations. Yet, once the quasicrystal is represented as a projection of a 
higher-dimensional Bravais lattice [2^, the phason flips are obtained from simple lattice translations along the extra 
dimensions. An important distinction between the new symmetry and phasons is that the generated transformation 
is continuous rather than discrete. Thus, there are no energy barriers between the degenerate shapes. Boundary 
conditions are expected to remove the degeneracy and will usually favor the symmetric and/or antisymmetric states 
[9|. The corresponding energy differences, however, should be exponentially small in nL. Therefore, we anticipate 
important consequences of the continuous degeneracy for the dynamics of ultra-thin fluid-supported films. 

On the face of it, the /C symmetry arises only from the translational invariance of the two-dimensional host system 
in which it is embedded. There is no apparent need for integrability of the host system. Yet, when one attempts 
to determine K. in practice for a non-integrable host system, the procedure of Eqs. p2)) - p6|) fails. The differential 
conditions imposed by these equations do not integrate to a well-defined function of {(f), h,p^,ph). The requirement 
that these conditions must define a definite K, function may then be a way of formulating the condition for integrability 
of the host system. 



The floating elastic sheets discussed here occur widely in physical systems, especially in molecular and nanoparticle 
monolayers 0, [l3,[30|. Thus we expect the degenerate shapes described above to have physical implications. First, 
one should in principle see the degeneracy in the model experiments on macroscopic sheets where this folding was 
initially demonstrated [Oj. However, only the symmetric and antisymmetric folds were reported. These sheets had 
lengths L not much greater than the width of the deformed region, so that boundary effects not considered here could 
have broken the degeneracy. In molecular sheets [3l| one may imagine a distortion of the y-invariant fold in which 
the asymmetry parameter varies slowly along the y direction, as demonstrated in Fig. [31 The extra energy for such 
a distortion must vanish as its wavelength goes to infinity, even though its amplitude is large. Here too, extra effects 
may act to select one or another of the asymmetric folds. For example, self- attraction favors the symmetric fold, since 
this variant approaches itself more closely than the others. On the other hand, the formation of a fold requires viscous 
dissipation, enhanced by shear stress. This stress is likely greater in the symmetric fold, so that the antisymmetric 
fold should be kinetically favored. 



,/ 







FIG. 3. (color online). Shape of a folded sheet whose y-invariance is broken by a "soft u-mode". Each xz cut is one of the 
degenerate solutions given by Eq. (|10|) with C2 = Q and c\ that changes linearly with y. The pressure is P = 0.5. 



The symmetry and degeneracy found here were an indirect and mysterious consequence of the integrability of the 
two-dimensional host system. Though we have discussed it in the context of a particular embedded system of physical 
interest, we expect that analogous symmetries occur elsewhere in the sG-mKdV hierarchy. These symmetries may 
ultimately deepen our understanding of integrability in these systems. In physical terms these degenerate deformations 
and the aforementioned "soft modes" associated with them may ultimately shed light on the experimental puzzles of 
how molecular folding and unfolding occur in practice |30l432l |. 



ACKNOWLEDGMENTS 

We thank Ron Lifshitz for a helpful discussion. This work was supported in part by the University of Chicago 
MRSEC program of the NSF under Award Number DMR 0820054. 



[1] B. Davidovitch, R. D. SchroU, D. Vella, M. Adda-Bedia, and E. A. Cerda, Proc. Natl. Acad. Sci. USA 108, 18227 (2011). 



[2] B. Davidovitch, R. D. SchroU, and E. Cerda, Phys. Rev. E 85, 066115 (2012). 

[3] H. King, R. D. SchroU, B. Davidovitch, and N. Menon, Proc. Natl. Acad. Sci. USA 109, 9716 (2012). 

[4] T. A. Witten, Rev. Mod. Phys. 79, 643 (2007). 

[5] S. T. Milner, J.-F. Joanny, and P. Pincus, Europhys. Lett. 9, 495 (1989). 

[6] E. Cerda and L. Mahadevan, Phys. Rev. Lett. 90, 074302 (2003). 

[7] D. VeUa, P. Aussillous, and L. Mahadevan, Europhys. Lett. 68, 212 (2004). 

[8] J. Huang, M. Juszkiewicz, W. H. de Jeu, E. Cerda, T. Emrick, N. Menon, and T. P. RusseU, Science 317, 650 (2007). 

[9] L. Pocivavsek, R. Dellsy, A. Kern, S. Johnson, B. Lin, K. Y. C. Lee, and E. Cerda, Science 320, 912 (2008). 
[10] B. D. Leahy, L. Pocivavsek, M. Meron, K. L. Lam, D. Salas, P. J. Viccaro, K. Y. C. Lee, and B. Lin, Phys. Rev. Lett. 

105, 058301 (2010). 
[11] J. Huang, B. Davidovitch, C. D. Santangelo, T. P. Russell, and N. Menon, Phys. Rev. Lett. 105, 038302 (2010). 
[12] D. P. Holmes and A. J. Crosby, Phys. Rev. Lett. 105, 038303 (2010). 
[13] H. Diamant and T. A. Witten, arXiv: 1009.2487] 
[14] B. Audoly, Phys. Rev. E 84, 011605 (2011). 

[15] H. Diamant and T. A. Witten, Phys. Rev. Lett. 107, 164302 (2011). 

[16] J. M. T. Thompson and G. W. Hunt, A General Theory of Elastic Stability (Wiley, London, 1973). 
[17] G. W. Hunt, M. K. Wadee, and N. Shiacolas, J. Appl. Mech. 60, 1033 (1993). 
[18] S. H. Lee and A. M. Waas, Int. J. Non-Linear Mech. 31, 313 (1996). 
[19] Q. Zhang and T. A. Witten, Phys. Rev. E 76, 041608 (2007). 
[20] B. Audoly and A. Boudaoud, J. Mech. Phys. Solids 56, 2401 (2008). 

[21] P. M. Reis, F. Corson, A. Boudaoud, and B. Roman, Phys. Rev. Lett. 103, 045501 (2009). 
[22] F. Brau, H. Vandeparre, A. Sabbah, C. Poulard, A. Boudaoud, and P. Damman, Nat. Phys. 7, 56 (2011). 
[23] H. Diamant and T. A. Witten, Talk presented at the APS Meeting, Boston, March 2012, 

http://meetings.aps.org/link/BAPS.2012.MAR.L52.10 
[24] M. Rivetti, C. R. Mec. 341, 333 (2013). 

[25] E. Cerda and L. Mahadevan, Proc. R. Soc. A 461, 671 (2005). 

[26] F. Gesztesy and H. Holden, Soliton Equations and Their Algebro-Geometric Solutions (Cambridge University Press, Cam- 
bridge, 2003), Chapter 2. 
[27] D. Levine, T. C. Lubensky, S. Ostlund, S. Ramaswamy, P. J. Steinhardt, and J. Toner, Phys. Rev. Lett. 54, 1520 (1985). 
[28] P. A. Kalugin, A. Yu. Kitaev, and L. S. Levitov, JETP Lett. 41, 145 (1985). 
[29] R. Lifshitz, Isr. J. Chem. 51, 1156 (2011). 
[30] K. Y. C. Lee, Ann. Rev. Phys. Chem. 59, 771 (2008). 
[31] Y. Zhang and T. M. Fischer, J. Phys. Chem. B 109, 3442 (2005). 
[32] A. Copal, V. Belyi, H. Diamant, T. A. Witten, K. Y. C. Lee, J. Phys. Chem. B 110, 10220 (2006). 



